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To  fabricate  micro-thermoelectric  generators  (pTEGs),  one  must  design  the  optimal  structure  of  the 
pTEGs  and  achieve  thermoelectric  thin  films  with  excellent  properties.  This  study  investigated  the  role 
of  the  dimensions  of  pTEGs,  including  the  length  of  the  thermoelements,  thickness  of  the  substrates, 
and  cross-sectional  area  of  the  thermoelements.  To  evaluate  the  power  generated  by  pTEGs  and  their  effi¬ 
ciency,  three-dimensional  models  of  pTEGs  were  subjected  to  finite  element  analysis.  Three-dimensional 
models  are  more  accurate  than  one-dimensional  models,  since  the  directions  of  the  heat  flux  and  electri¬ 
cal  current  are  not  parallel  in  pTEGs.  The  governing  equations  were  derived  from  the  Seebeck  effect  and 
Peltier  effect,  which  show  thermoelectric  energy  conversion.  In  the  simulation,  the  substrate,  n-type 
material,  and  p-type  material  were  assumed  to  be  silicon,  Bi2Te3,  and  Sb2Te3,  respectively.  We  calculated 
the  thermoelectric  power  generated  by  the  pTEGs  and  their  thermoelectric  energy  conversion  efficiency. 
These  two  evaluation  indices  represent  the  performance  of  pTEGs.  The  thermoelectric  simulation  pro¬ 
duced  design  guidelines  for  high-performance  pTEGs. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

For  the  sustainable  development  of  humankind  and  to  stop  cli¬ 
mate  change  as  outlined  in  the  Kyoto  Protocol,  the  use  of  conven¬ 
tional  energy  sources  such  as  fossil  fuels  will  have  to  be  limited  in 
the  very  near  future.  New  alternative  sources  of  energy  such  as 
wind,  hydropower,  and  solar  energy  are  important  research  topics 
in  physics,  chemistry,  biology,  and  mechanics.  As  one  such  alterna¬ 
tive,  thermoelectric  energy  has  many  favorable  characteristics: 
light  weight,  small  scale,  and  low  manufacturing  cost  [1].  Given 
these  advantages,  thermoelectric  modules  have  numerous  applica¬ 
tions,  such  as  in  the  energy  source  for  autonomous  microsystems, 
harvesting  energy  from  the  waste  heat  from  high-performance 
microprocessors,  body-powered  wearable  and  pocket  devices, 
and  semipermanent  power  for  implanted  medical  devices  [2-5]. 

Thermoelectric  effects  involve  the  reversible  conversion  of 
energy  from  thermal  energy  to  electrical  energy.  The  Seebeck  effect 
refers  to  forward  energy  conversion  and  the  Peltier  effect  is  reverse 
conversion.  These  effects,  which  were  discovered  about  200  years 
ago,  have  not  been  used  in  power  generator  due  to  the  low  energy 
conversion  efficiency  [1].  However,  many  researchers  have  studied 
thermoelectric  power  generation  since  the  possibility  of  enhancing 
the  thermoelectric  conversion  efficiency  using  nanotechnology 
was  reported  [6,7]. 

To  fabricate  high-performance  micro-thermoelectric  generators 
(pTEGs),  thermoelectric  materials  deposited  using  microelectro¬ 
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mechanical  systems  (MEMS)  processes  in  the  form  of  thin  films 
should  have  excellent  thermoelectric  properties,  and  a  great  deal 
of  active  research  is  now  examining  high-performance  materials 
[8,9].  In  addition  to  high-quality  thermoelectric  materials,  the  opti¬ 
mal  design  of  an  integrated  module  that  includes  thermoelectric 
materials,  electrodes,  and  substrates  is  indispensable  to  achieve 
high-performance  pTEGs.  Some  researchers  have  already  studied 
the  geometric  dependence  of  thermoelectric  modules  using  one¬ 
dimensional  models  [10-13].  However,  these  studies  cannot 
explain  the  effects  of  various  geometric  shapes.  Moreover,  an  accu¬ 
rate  analysis  is  impossible  when  the  directions  of  the  heat  flux  and 
electrical  current  are  not  parallel  in  the  thermoelements  of  pTEGs 
made  from  thin  thermoelectric  materials  compared  to  the  cross- 
sectional  area.  In  this  study,  we  investigated  the  optimal  structure 
of  high-performance  pTEGs  using  the  finite  element  method  with 
three-dimensional  models. 

2.  Formulation 

2A.  Geometric  configuration 

We  considered  a  conventional  thermoelectric  unit  composed  of 
a  single  p-n  pair.  Fig.  1  shows  the  7t-shaped  thermoelectric  unit  in 
which  the  p-  and  n-type  thermoelements  are  connected  in  parallel 
thermally  and  in  series  electrically.  The  electrodes  of  the  p-  and  n- 
type  thermoelements  are  connected  electrically.  In  addition,  the 
upper  and  lower  substrates  play  a  role  in  thermal  conduction  from 
heat  sources  and  cooling  parts  to  the  electrodes  linking  the  ther¬ 
moelements.  They  also  function  as  an  electrical  insulation  layer 
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between  the  outer  heat  source  and  cooling  parts.  We  assume  that 
the  pTEGs  are  constructed  from  such  thermoelectric  units  joined  in 
series. 

2.2.  Governing  equations 

From  the  finite  element  analysis,  we  obtain  steady-state  solu¬ 
tions  for  the  governing  equations  explained  below.  The  governing 
equations  express  the  thermal,  electrical,  and  thermoelectric  ef¬ 
fects  in  the  thermoelectric  materials. 

When  temperature  gradient  VT  is  given,  the  electrical  field  E 
generated  in  the  thermoelectric  material  is  given  by 

E  =  ocVT-pl  (1) 

from  the  Seebeck  effect  and  Ohm’s  law.  Here,  a,  p,  and  J  denote  the 
Seebeck  coefficient,  resistivity,  and  current  density,  respectively. 
The  Peltier  effect,  which  explains  the  heat  flux  q  in  the  thermoele¬ 
ments,  is  expressed  by 

q  =  ocTJ  +  tcVT.  (2) 

Here,  the  second  term  explains  the  thermal  conduction  with 
thermal  conductivity  k.  Eqs.  (1)  and  (2)  are  partial  differential 
equations  that  explain  the  thermoelectric  effects  in  thermoelectric 
materials 

Vq=JE,  (3) 

V-J=  0,  (4) 

E=-VV,  (5) 

express  the  energy  conservation  law,  continuity  of  current,  and  def¬ 
inition  of  the  electrical  potential  V,  respectively.  Since  we  have  five 
equations  and  five  unknown  variables,  we  can  obtain  the  solution 
by  solving  these  simultaneous  partial  differential  equations  with 
the  finite  element  method. 

Note  that  Eqs.  (1)  and  (2)  show  Ohm’s  law  and  Fourier’s  law 
when  a  =  0.  In  the  electrodes,  since  the  governing  equations  consist 
only  of  electrical  and  thermal  conduction  equations,  Eqs.  (1)  and 
(2)  are  used  in  the  analysis.  Consequently,  we  can  simply  calculate 
the  temperature  distribution  in  the  substrates  using  Eq.  (2). 

2.3.  Simulation 

We  assumed  the  use  of  p-  and  n-type  materials,  such  as  Sb2Te3 
and  Bi2Te3,  which  have  excellent  thermoelectric  properties  at  room 


Table  1 

Properties  of  the  materials  used  in  the  finite  element  analysis. 


Materials 

Thermal 

conductivity 

(K) 

Electrical 

conductivity 

(P) 

Seebeck 
coefficient  (a) 

n-Type  thermoelement 

1.6  W/m  I< 

7.69  x  104  S/m 

-2.28  x  10“4V/I< 

(Bi2Te3) 

p-Type  thermoelement 

2.1  W/m  K 

9.62  x  104  S/m 

1.71  x  10“4V/K 

(Sb2Te3) 

Electrode  (Cu) 

Substrate  (Si) 

350  W/m  K 
130  W/m  K 

5.90  x  108  S/m 

temperature.  The  substrates  and  electrodes  are  silicon  and  copper, 
respectively,  which  are  widely  used  materials  in  the  MEMS  pro¬ 
cess.  The  properties  of  these  materials  are  shown  in  Table  1  [13]. 

The  boundary  conditions  are  described  in  Fig.  1.  The  current 
flows  through  one  side  in  an  electrode,  and  the  electrical  potential 
is  zero  at  the  surface  of  the  counter-side  electrode.  The  tempera¬ 
tures  at  the  top  of  the  upper  substrate  and  bottom  of  the  lower 
substrate  are  300  K  and  285  K,  respectively.  For  application  as  har¬ 
vesters  of  wasted  heat  for  wearable  and  pocket  electronic  devices, 
we  performed  the  analysis  with  a  relatively  small  difference  in 
temperature  [4].  The  other  boundaries  are  assumed  to  be  adiabatic 
boundary  conditions. 

For  a  given  electrical  current  flow  into  the  unit,  we  calculated 
the  electrical  potential,  temperature  distribution,  and  heat  flux. 
Then,  the  efficiency  r\  of  this  unit  is  determined  as 

ri  =  P/d  in,  (6) 

where  P  is  the  power  generated  in  the  unit  and  Qm  is  the  amount  of 
heat  flux  through  the  hot  side  surface  of  the  unit. 


3.  Results  and  discussion 

First,  we  investigated  the  effect  of  the  substrate,  as  shown  in 
Figs.  2  and  3.  As  the  thickness  of  the  substrate  increases,  the 
thermal  loss  in  the  substrate  becomes  larger  than  that  in  the  ther¬ 
moelements.  For  this  reason,  the  maximum  difference  in  tempera¬ 
ture  and  the  maximum  electrical  potential  of  the  unit  are  small 
when  the  substrates  are  thick.  Therefore,  the  generated  power  de¬ 
creases  as  the  substrate  gets  thicker.  Although  the  heat  flux 
through  the  hot  side  surface  also  decreases  in  thick  substrates,  it 


Fig.  2.  Maximum  power  generated  in  thermoelectric  units  with  various  substrate 
thicknesses. 
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Fig.  3.  Maximum  efficiency  obtained  in  thermoelectric  units  with  various  substrate 
thicknesses. 


is  small  compared  to  the  reduction  in  the  power  generated  by  thick 
substrates.  Hence,  a  unit  with  thick  substrates  has  low  efficiency. 

To  determine  the  geometric  effects  of  the  thermoelements,  the 
dependences  of  the  power  and  efficiency  on  the  length  of  the  ther¬ 
moelements  are  shown  in  Figs.  4  and  5,  respectively.  Generally, 
power  generated  in  the  unit  increases  when  the  thermoelements 
are  small  due  to  the  decrease  in  the  electrical  resistance.  However, 
for  short  thermoelements,  the  difference  in  temperature  between 
the  upper  and  lower  sides  of  the  thermoelements  decreases.  For 
this  reason,  the  length  of  the  thermoelements  is  optimized  to  give 
the  highest  power.  The  efficiency  tends  to  increase  with  the  length 
of  the  thermoelements  due  to  the  dominant  effect  of  the  heat  flux 
decrement. 

Finally,  we  investigated  the  characteristics  of  a  unit  with  ther¬ 
moelements  having  various  cross-sectional  areas.  The  open  circuit 
voltages,  which  are  directly  related  to  the  maximum  difference  in 
temperature  from  Eq.  (1),  are  large  for  units  with  small  cross-sec¬ 
tional  areas.  However,  the  short  circuit  current  is  small,  since  the 
electrical  resistance  is  large  when  the  cross-sectional  area  of  the 
thermoelements  is  small.  Consequently,  a  large  power  output 
is  obtained  in  a  unit  with  thermoelements  having  large  cross- 
sectional  areas,  as  shown  in  Fig.  6.  However,  the  efficiency  de¬ 
creases  as  illustrated  in  Fig.  7,  when  the  cross-sectional  area  is 


Fig.  4.  Maximum  power  generated  in  the  thermoelectric  units  with  thermoele¬ 
ments  of  various  lengths. 


Fig.  5.  Maximum  efficiency  obtained  in  thermoelectric  units  with  thermoelements 
of  various  lengths. 


Fig.  6.  Maximum  power  generated  in  thermoelectric  units  with  thermoelements  of 
various  cross-sectional  areas  (Lx  L  pm2). 


Fig.  7.  Maximum  efficiency  obtained  in  thermoelectric  units  with  thermoelements 
of  various  cross-sectional  areas  (L  x  L  pm2). 
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Fig.  8.  Maximum  power  per  unit  area  obtained  in  thermoelectric  units  with 
thermoelements  of  various  cross-sectional  areas  (L  x  L  pm2). 


large  because  the  heat  flux  through  the  hot  side  is  proportional  to 
the  cross-sectional  area  of  the  thermoelements.  For  the  same  rea¬ 
son,  the  power  generated  per  unit  area  is  small  in  a  unit  with  a 
small  cross-sectional  area,  as  shown  in  Fig.  8.  In  contrast,  the 
dependence  of  the  power  output  per  unit  area  is  similar  to  that 
of  the  generated  power  on  the  substrate  thickness  and  thermoele¬ 
ment  length  for  the  same  unit  area.  The  method  used  to  evaluate 
the  performance  of  pTEGs  can  be  summarized  using  the  generated 
power,  power  output  per  unit  area,  and  efficiency.  The  efficiency 
given  by  Eq.  (6)  is  the  ideal  parameter.  However,  measuring  the 
efficiency  of  a  thermoelectric  module  in  practice  is  difficult  due 
to  the  difficulty  in  measuring  the  heat  flux.  For  this  reason,  pTEG 
manufacturers  provide  the  generated  power  or  power  output  per 
unit  area  in  their  data  sheets.  To  develop  pTEGs  considering  the 
usage  and  heat  source,  estimating  the  efficiency  from  an  analysis 
of  the  module  is  essential. 


4.  Conclusions 

The  dependence  of  thermoelectric  performance  on  the  sub¬ 
strate  thickness  and  the  length  and  cross-sectional  area  of  the  ther¬ 
moelements  was  investigated  using  finite  element  analysis  with  a 
three-dimensional  model.  As  the  substrate  gets  thicker,  the  ther¬ 
moelectric  performance  deteriorates  due  to  thermal  loss  from  the 
substrate.  The  thermoelements  have  an  optimal  length  with  the 
highest  power.  High  efficiency  is  obtained  when  the  length  of  the 
thermoelements  is  large.  The  power  generated  declines  with  the 
cross-sectional  area  of  the  thermoelements,  while  the  efficiency 
shows  the  opposite  trend.  The  exact  prediction  of  the  efficiency 
using  the  analysis  helps  the  engineer  evaluate  the  thermoelectric 
module.  To  design  high-performance  pTEGs,  the  design  parameters 
should  be  optimized  by  analyzing  a  three-dimensional  model. 
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